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Abstract. We present an analysis of the spatial fluctua- 
tions of the 2-10 keV Cosmic X-ray Background (CXB) as 
measured from 22 high galactic latitude (|b| > 25°) fields 
observed with the MECS instrument on-board BeppoSAX. 
This technique allowed us to probe extragalactic source 
counts a factor 3-4 fainter than is possible with direct 
measurements of pointlike sources in MECS deep fields. 
The slope of the 2-10 keV logA'-logS' relationship is found 
to be still close to the "Euclidean" one (7 = 1.5) down 
to our flux limit of ~ 1.5 x 10"^** erg cm^-^ s ^ , where 
the contribution of discrete sources to the 2-10 keV CXB 
amounts to ~ 40 — 50%. Source counts derived from the 
analysis presented in this letter are in very good agree- 
ment both with those directly measured with ASCA and 
BeppoSAX deep surveys at bright fluxes and with a first 
estimation of the 2-10 keV Chandra logA^-log5' at fainter 
fluxes. 
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1. Introduction 

The understanding of the origin of the Cosmic X-ray Back- 
ground (CXB, Giacconi et al. 1962) has been and still is 
one of the main topics of X-ray astronomy. It is now widely 
accepted that the CXB is the results of the superposition 
of the emission of faint extragalactic discrete sources and 
many surveys of the X-ray sky have been carried out to 
directly resolve the CXB and study the nature of these 
sources. 

The soft X-ray energy band (0.5-2 keV) has been ex- 
tensively investigated with ROSAT with which a large 
fraction of the soft CXB (70 — 80%) has been resolved into 
discrete sources (Hasinger et al. 1993, 1998). Spectroscopic 
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identifications of their optical counterparts have revealed 
that the large majority (> 80%) are AGNs, mostly with 
broad emission lines (Schmidt et al. 1998). 

The harder 2-10 keV sky has become accessible to deep 
surveys only over the past few years thanks to the imaging 
instruments of ASCA and BeppoSAX. About 30% of the 
CXB in this band has been directly resolved into discrete 
sources and the logiV-log5 relationship has been derived 
down to a flux limit of 5 x 10^^^ erg cm~^ s ^ (Ueda et 
al. 1998; Cagnoni, Delia Ceca & Maccacaro 1998; Giommi 
et al. 1998; Della Ceca et al. 1999; Ueda et al. 1999a, 
1999b; Fiore et al. 1999; Giommi, Fiore & Perri 1999; Fiore 
et al. 2000a; Giommi, Perri & Fiore 2000). 

A number of deeper surveys carried out with Chandra 
are now starting to probe the 2-10 keV CXB at much 
fainter fluxes. However, at least for some time, these can 
only cover a small area of sky with the consequence that 
the statistics on source counts (and optical identifications) 
is still very limited (Brandt et al. 2000; Fiore et al. 2000b; 
Mushotzky et al. 2000; Hornschemeier et al. 2000). 

In this letter we present the results of a fluctuation 
analysis of BeppoSAX MECS images performed in the 2- 
10 keV energy band that allowed us to study the logiV- 
logS* down to fluxes of ~ 1.5 x lO"""^"* erg cm~^ s . The 
analysis of the spatial fluctuations of the CXB is a pow- 
erful method that allows us to investigate the logA^-logS' 
relationship beyond the flux limit where single sources be- 
came too faint, or confused, to be directly detected. The 
theory behind this technique has been described for the 
first time by radio astronomers (Scheuer 1957, 1974; Con- 
don 1974) and many applications of this method to X- 
ray data have successfully predicted source counts at faint 
fluxes (Hamilton & Helfand 1987; Barcons & Fabian 1990; 
Hasinger et al. 1993; Barcons et al. 1994; Butcher et al. 
1997; Gendreau, Barcons & Fabian 1998). Basically, the 
method consists of i) comparing the observed spatial fluc- 
tuations distribution of the CXB with a set of model dis- 
tribution curves corresponding to different hypotheses on 
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the power law slope 7 and normalization K of the logA''- 

log5 relationship and ii) by means of a fitting procedure 
deriving best fit parameters on source counts. 

The main contribution to the observed CXB spatial 
fluctuations comes from sources at a Q\ix level S corre- 
sponding to a source density of about 1-2 sources per beam 
(Scheuer 1974). This condition determines the sensitivity 
limit of this technique. 

2. Data analysis 

The instruments aboard the BeppoSAX Satellite (Boella 
et al. 1997a) include a Low Energy (0.1-10 keV) Con- 
centrator Spectrometer (LEGS, Parmar ct al. 1997), a 
Medium Energy (1.3-10 keV) Concentrator Spectrometer 
composed by three units (MECS, Boella et al. 1997b), a 
High Pressure Gas Scintillation Proportional Goimtcr (3- 
120 keV) (HPGSPG, Manzo et al. 1997) and the Phoswich 
Detector System (15-300 keV) (PDS, Prontera et al. 1997). 
All these instruments point in the same direction and arc 
collectively called Narrow Field Instruments (NFIs). In 
addition BeppoSAX carries two wide field (20 x 20 de- 
grees FWHM) instruments operating in the 2-30 keV en- 
ergy band (Wide Field Cameras, WFC, Jager et al. 1997) 
and pointing in diametrically opposed directions perpen- 
dicular to the NFIs. 

The data used in our fluctuation analysis consist of 
22 non-overlapping high galactic latitude (|b| > 25°) 
MECS fields pointed at "blank" parts of the sky, that 
is not centered on previously known X-ray targets. About 
half (13) of these fields are Secondary Pointing NFIs ob- 
servations, produced when the primary target (usually 
the Galactic Center) was observed with one of the Wide 
Field Cameras. The other observations include 3 non- 
overlapping fields centered near Polaris, the BeppoSAX 
default pointing position in case of safe mode, 2 follow- 
up observations of Gamma Ray Bursts (GRBs) carried 
out early in the mission in which X-ray afterglows of the 
GRBs were not detected, 3 images, taken from the pub- 
lic archive, not centered on known X-ray sources, and 1 
image taken during the Leonids meteorites crossing on 
November 17, 1999, when the satellite was oriented so as 
to minimize the probability of impact with one of the me- 
teorites. The average cxposmc of the 22 fields is ^ 88 ks 
while typical values of galactic Nh, as derived by the 
21 cm measurements of Dickey & Lockman (1990), are 
~ 4 X 10^" cm^^. No point sources with 2-10 keV flux 
greater than 5 x 10^^^ erg cm~^ s ^ are detected in any 
of these fields. 

As discussed above, the MECS instruments is com- 
posed of three nearly identical units. One of these 
(MECSl) failed on May 7th 1997. To ensure uniformity we 
have used data only from the imits MECS2 and MECS3, 
to which we refer hereafter as the MECS instrument. To 
avoid difficulties due to complex variations in sensitivity 
across the field-of-view (FOV) we have restricted our anal- 




Fig. 1. MECS 2-10 keV image of one of the "blank" fields used 
in our analysis. The 4 equal quadrants in which the central 8 
arcminutes circular region has been divided are also shown. 



ysis to the central 8 arcminutes circular regions (see Figure 
1) where the sensitivity is higher and vignetting and Point 
Spread Function (PSF) variations are small. 

In order to enhance the sensitivity of our analysis as 
well the statistics of the CXB granularity each circular re- 
gion has been divided in 4 equal quadrants (each covering 
an area of ~ 50.3 arcmin^, see Figure 1) for a total of 88 
independent measurements of the CXB. This sets our sen- 
sitivity to 5 ~ 1.5 X 10~^^ erg cm~^ s ^, the exact value 
depending on the details of the logA^'-logS' relationship. At 
this flux level one source corresponds to ^ 10 — 15 MECS 
net counts, well above the "1 photon per source limit" 
which is considered as a further limit of the fluctuation 
analysis method. 

2.1. The Non X-ray Background 

To estimate the amount of cosmic X-rays we must subtract 
the non X-ray background (NXB, that is the instrument 
noise plus the particle background) which, in the central 
8 arcminutes of the MECS, accounts for about half of the 
total counts. 

To measure the NXB we have used, as in Parmar et al. 
(1999) and in Vecchi et al. (1999), BeppoSAX MECS "dark 
Earth fields" , that is MECS images taken when the point- 
ing direction is occulted by the part of the Earth that is 
not illuminated by the Sun. We have accumulated "dark 
Earth" observations from July 1996 through September 
1999 for a total exposure time of 3.85 Ms. Counts between 
channels 44 and 200 (the energy range used to estimate 
the CXB, see below) have been extracted from the inner 
8 arcminutes circular regions of the images. The cor re- 
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Fig. 2. MEGS Non X-ray Background count rate (channels 
44-200, central 8 arcminutes circular regions) as a function of 
time. "Dark Earth" data have been accumulated from 1996 
July 6 (t = 0) to 1999 September 11 for a total exposure time 
of 3.85 Ms. 

spending MECS light curve is shown in Figure 2. A grad- 
ual decline of the count rate with time is evident. From 
July 1996 to September 1999 the amount of this reduction 
is ~ 15 — 20%. We modeled this behaviour with the law 
cr{t) ~ cr(0) — 2.07 x 10^^^ x t, where t is time in seconds 
and cr(0) ~ 1.11 x 10~^ s~^ is the NXB count rate value 
at the beginning of the mission (t = 0, 1996 July 6). A 
similar instrumental background reduction has been re- 
ported for the LEGS instrument on-board BeppoSAX by 
Parmar et al. (1999). 

We also investigated potential systematic errors in the 
estimation of the NXB, such as significant deviations of 
count rate from the law given above or the statistical pois- 
son noise of the NXB counts in a MECS quadrant. In both 
cases, given the typical exposure times of our MECS fields 
(T ~ 90 ks), such effects are negligible with respects to the 
measured fluctuations of the net counts. 

2.2. The CXB flux distribution 

Net counts between channels 44 and 200 (2.0-9.0 keV) 
have been extracted from the central 8 arcminutes quad- 
rants of the 22 "blank" fields and converted to fluxes 
in the 2-10 keV band. A power law spectrum with en- 
ergy index a = 0.6 has been assumed, as indicated by 
ASCA and BeppoSAX deep 2-10 keV X-ray surveys (Ueda 
et al. 1999b; Giommi, Perri & Fiore 2000). Corrections 
for the Galactic absorbing neutral hydrogen column den- 
sity Nh have also been applied. Counts in the channel 
range 200-220 (9.0-10.0 keV) have been excluded since the 
NXB dominates this energy band. The flux distribution 
of the 88 independent measurements of the Cosmic X- 
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Fig. 3. Distribution of the 2-10 keV Gosmic X-ray Background 
measurements (histogram) . The best fit model curve from our 
ML estimation (solid line) is also shown. 

ray Background is shown in Figure 3. The average flux is 
~ 2.3 X 10"" erg cm^^ s"Vleg"^, a value that is - 35% 
higher than that of HEAOl (Marshall et al. 1980) and 
slightly higher (~ 15%) than the ASCA GIS value (Miyaji 
et al. 1998) but in good agreement with the Wisconsin 
rocket data (McCammon & Sanders 1990 and references 
therein) and with the MECS measurement obtained from 
a smaller data sample by Vecchi et al. (1999). 

3. Fluctuation Analysis 

To compute the predicted CXB spatial fluctuations we 
have adopted an analytical approach that, besides the 
statistical fluctuations due to the expected number of X- 
ray sources, includes coimting statistics and the effects of 
MECS vignetting and PSF. 

Our method can be summarized as follows. A power 
law form for the logA^-log5 relationship is assumed, N{> 
S) = K (S/So)^^ , where K is the number of sources per 
unit area with flux S > So- 

The expected number TV of sources in one beam (one 
MECS quadrant) in the flux range from S to S + dS is 
then obtained simply multiplying the differential of the 
logTV-logS' with the area of a MECS quadrant. Next, the 
probability of a fluctuation of the number of sources is 
computed assuming Poisson statistics and the correspond- 
ing photon counts distribution. Counting statistics and vi- 
gnetting effects are included in the procedure by means of 
two independent probability convolutions with i) the Pois- 
son distribution of counts themselves and ii) the off-axis 
angles probability distribution of sources in the MECS 
FOV. Probability convolutions have been computed by 
means of a Fast Fourier Transform (FFT) technique. 
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Fig. 4. 2-10 keV integral source counts. Solid and dashed "bow tie" regions are 90% and 68% contours from the present 

BeppoSAX fluctuation analysis. Dotted lines are the 1 a contours from the ASCA fluctuation analysis of Gendreau, Barcons & 
Fabian (1998). Filled circles represent the -BeppoSAX 2-10 keV logA^-logS relationship (Giommi, Perri & Fiore 2000) and open 
circles are Chandra source counts from the SSA13 field (Mushotzky et al. 2000). 



The 2-10 kcV MECS PSF size (the radius where 80% 
of the photons are collected is ~ 2.6 arcminutes, Boella 
et al. 1997b) is non-negligible compared to the regions 
were we extract counts and consequently photons due to 
sources just outside the 8 arcminutes quadrant cannot be 
neglected. Similarly to the ASCA fluctuation analysis of 
Gendreau, Barcons & Fabian (1998) we have then divided 
a MECS quadrant and "nearby regions" with 2 arcmin- 
utes size boxes and computed the fraction of photons that 
fall inside a quadrant at different FOV positions. We have 
next constructed the off-axis angles probability distribu- 
tion of sources and included this effect with a further con- 
volution. This procedure is repeated for all fluxes between 
S'min, the flux at which the logiV-logS* integrated inten- 
sity exceeds the CXB value (as estimated in Vecchi et al. 
1999), and S'max = 5 x lO^"'^'' erg cm^^ s ^ (sources with 
S > S'max are absent from our data and do not contribute 
to the CXB fluctuations). Finally, from the probability 
convolution of the computed flux distributions wc obtain 
the expected distribution of CXB spatial fluctuations cor- 
responding to the assumed values of 7 and K. 

We have compared the observed CXB spatial fluc- 
tuations with a set of model curves computed assum- 
ing different trial logA''-log5 relationships. A maximum 



likelihood (ML) test has been used to estimate the best 
flt, 68% and 90%, {AS = 2.3 and 4.6) constraints on 
the 2-10 keV logA^'-logS' power law slope 7 and normal- 
ization K (as defined by equation (2) with 6*0 = 1 x 
10~^^ erg cm~^ s ^). The best fit values are 7 = 1.5 ± 0.2 
and K = SSS^gg^ deg~^ where the quoted errors corre- 
spond to 90% confidence level for one interesting parame- 
ter {AS = 2.7). The solid line in Figure 3 shows the model 
best-fit fiux distribution which approximates well the ob- 
served one (a Kolmogorov-Smirnov test gives a probability 
> 95% that the two distributions are equal). 

The 90% and 68% constraints of our fluctuation analy- 
sis on the 2-10 keV integral logA^-logS* arc shown in Figure 
4. For comparison the 1 a contours from the ASCA SIS 
fluctuation analysis of Gendreau, Barcons & Fabian (1998) 
and 2-10 keV source counts from BeppoSKK. (Giommi, 
Perri & Fiore 2000) and Chandra (Mushotzky et al. 2000) 
surveys are also plotted. 

4. Simulations 

In order to check our procedure we have carried out exten- 
sive simulations of MECS images using the data simulator 
available at the SeppoSAX Science Data Center (Giommi 
& Fiore 1997). This tool fully includes MECS instrumen- 
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tal features as telescope vignetting, PSF shape, non X-ray 
background and Poisson noise. 

We generated 3 sets of one hundred MECS images with 
pointlike sources following 3 different logA/'-logS' distri- 
butions. Sources have been simulated in the flux range 
'S'min — Smax (sce above) with exposure time T = 90 ks. 

Each set of simulated images have been analyzed fol- 
lowing the same procedure used for the real data. The 
resulting flux distributions have then been compared with 
the model curves and the maximum likelihood test per- 
formed. 

We have found that the input values for the power law 
slope 7 and the normalization K of the logA''-log5 are in 
all 3 cases within the 68% confidence contours of the fitted 
ones indicating that no significant bias in their estimation 
affects our analysis. 

5. Summary and Conclusions 

We have performed a CXB fluctuation analysis based 
on 22 non-overlapping high galactic latitude (|b| > 25°) 
MECS images of "blank" parts of the sky. The total 
exposure time is ~ 1.94 Ms while the area covered is 
~ 1.23 deg^. The analysis of the CXB spatial fluctuations 
has allowed us to extend the 2-10 keV logiV-logS* rela- 
tionship down to fluxes ~ 1.5 x 10^^^ erg cm~^ s ^. At 
this flux level the contribution of discrete sources to the 
2-10 keV Cosmic X-ray Background is - 40 - 50%. The 
best fit values for the slope 7 and the normalization K of 
the logA^-logS' are 7 = 1.5 ± 0.2 and K = 336^^^'' dcg"^ 
(quoted errors correspond to 90% confidence level). Fig- 
ure 4 summarizes our results. A good agreement with the 
ASCA fluctuation analysis of Gendreau, Barcons & Fabian 
(1998) and with BeppoSAX direct source counts at bright 
fluxes (Giommi, Perri & Fiore 2000) and a first estimate of 
the Chandra 2-10 keV logAf-logS (Mushotzky et al. 2000) 
at faint fluxes is found. 
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